Revised 9 February 2004 
FTUAM 03-19 
hep-ph/0312187 



Regge analysis of pion-pion (and pion-kaon) scattering 
for energy s 1/2 > 1.4 GeV 



J. R. Pelaez 

Departamento de Fisica Teorica, II 
(Metodos Matemdticos) , 
Facultad de Ciencias Fisicas, 
Universidad Complutense de Madrid, 
E-28040, Madrid, Spain. 

and 

F. J. Yndurain 

Departamento de Fisica Teorica, C-XI 
Universidad Autonoma de Madrid, 
Canto Blanco, 
E-28049, Madrid, Spain. 



Abstract 

We perform a detailed Regge analysis of NN, irN, KN, tttt and irK scat- 
tering. From it, we find expressions that represent the nn scattering amplitudes 
with an accuracy of a few percent, for exchange of isospin zero, and ~ 15% for ex- 
change of isospin 1, and this for energies s 1//2 > 1.4 GeV and for momentum trans- 
fers [tl 1 / 2 < 0.4 GeV. These Regge formulas are perfectly compatible with the 
low energy (s 1//2 ~ 1.4 GeV) scattering amplitudes deduced from nn phase shift 
analyses as well as with higher energy (s 1 ^ 2 > 1.4 GeV) experimental nn cross 
sections. They are also compatible with NN, KN and ttN experimental cross 
sections using factorization, a property that we check with great precision. This 
contrasts with results from current phase shift analyses of the nn scattering am- 
plitude which bear little resemblance to reality in the region 1.4 < s 1 ^ 2 < 2 GeV, 
as they are not well defined and increasingly violate a number of physical re- 
quirements when the energy grows. nK scattering is also considered, and we 
present a Regge analysis for these processes valid for energies s 1 ^ 2 > 1.7 GeV. 

As a byproduct of our analysis we present also a fit of AW, nN and KN 
cross sections valid from cm. kinetic energy -Ekin — 1 GeV to multi TeV energies. 
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-REGGE ANALYSIS OF PION-PION (AND PION KAON) SCATTERING - 

1. Introduction 

A precise and reliable knowledge of the tttt scattering amplitude has become increasingly important in the last 
years. This is so, in particular, because tttt scattering is one of the few places where one has more observables 
than unknown constants in a chiral perturbation theory analysis, so it provides a window to higher order 
terms. Moreover, an accurate determination of the S wave scattering lengths, and of the phase shifts at 
s i/2 _ mK ^ provide essential information for two subjects under intensive experimental investigation at 
present, viz., pionic atom decays and CP violation in the kaonic system. In recent papers, Ananthanarayan, 
Colangclo, Gasser and LeutwylerW (that we will denote by ACGL), Colangclo, Gasscr and Leutwyler,! 2 ' 
Descotes et alJ 3 ' and Kamihski, Lesniak and Loiseau^ have used experimental information, analyticity and 
unitarity (in the form of the Roy equations) and, in ref. 2, chiral perturbation theory, to construct the mr 
scattering amplitude at low energy, s 1 / 2 < 0.8 GeV. For these analyses one needs as input the imaginary 
part of the mr amplitudes above the energy at which the Roy analysis stops; in particular, one needs the 
scattering amplitudes for s 1 / 2 above 1.4 GeV, which will be the subject of the present paper. 

Unfortunately, the authors in refs. 2, 3 take their tttt scattering amplitude in this energy region 
from ACGL,M which presents a number of serious drawbacks. 1 First of all, the input scattering amplitude 
at energy s 1 / 2 > 2 GeV which these authors use (following Pennington ^) is not physically acceptable, 
as it contradicts known properties of standard Regge theory and, moreover, is quite incompatible with 
experimental 2 tttt total cross sections, ^ and this in spite of the large errors assumed by ACGL. Secondly, 
the scattering amplitude for 1.4 GeV < s 1 / 2 < 1.9 GeV that ACGL (and, following them, the authors in 
refs. 2, 3) use is obtained from phase shift analyses, specifically the Cern-Munich set of analyses, ^ which 
are subject to large uncertainties and which, indeed, can be shown to contradict a number of physical 
requirements. [Although we will not discuss this here (see ref. 9), it is also clear that the errors ACGL, and 
the authors in ref. 2, take for some of their lower energy experimental input data are excessively optimistic 
and, moreover, certain of their chiral parameters are likely to be biased' 10 ']. One should imagine that the 
use of incorrect high energy input should lead to inconsistent low energy output. In fact, this occurs in the 
work by Colangelo, Gasser and Leutwyler,! 2 ! where the central values are probably displaced and the errors 
claimed are excessively optimistic and lead to several mismatches, as shown in refs. 9, 11. 

In the present note we will not concern ourselves with the reliability or otherwise of the low energy 
consequences of faulty high energy input, but will concentrate our efforts in ascertaining what a correct 
high energy input should be. To do this, we will perform a detailed Regge analysis and show that it is 
compatible with experimental data for all values of s 1 / 2 > 1.4 GeV (for some tttt processes, even down to 
s 1 / 2 ^ i GeV). The resulting tttt amplitudes, summarized in Eqs. (4), (5), (11) and (18) and Table II below, 
should provide a correct and accurate input for dispersive studies of tttt scattering. 

Our analysis will be an improvement on standard ones not only for tttt and ttK , but even for ttN, 
KN and NN in that we will be able to give an accurate description of the amplitudes for energies ranging 
from a kinetic enrgy in the center of mass E^ n ~ 1 GeV to the TeV region. This accuracy reaches the level 
of a very few percent for zero isospin exchange, and it is less precise for the isospin 1 exchange amplitude, 
for which the errors may go up to ~ 15% at low energy. 

An analysis of high energy ttK scattering is possible by a straightforward extension of the methods 
here; it is given in Sect. 3, where we present precise Regge formulas for zero isospin exchange, valid for 
energies s 1 / 2 > 1.7 GeV. 

The analysis of tttt, ttK scattering up to (relatively) low energies, <~ 14 GeV, is described in 
Sects. 2, 3; in Sect. 4 we extend it to multi TeV energies. As a byproduct of our analysis we present 

1 In ref. 4, the Regge parameters of ACGL are also used for ttK scattering; perhaps this is the reason why they are 
not able to get a satisfactory description of this process. 

2 It should be noted that Pennington has publicly stated (in the Conversano workshop, July 2003) that his analysis, 
tenable in 1974, is superseded by more recent developments, both experimental and theoretical. In fact, already 
by 1977 it was clear to experts that standard Regge behaviour also holds for tttt scattering; see, e.g., Froggatt and 
Petersen, ref. 6, who use the correct Regge behaviour in their dispersive analysis of tttt scattering. 
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also a parametrization of NN, irN and KN total cross sections compatible with the Froissart bound and 
valid from i?kin — 1 GeV to <~ 30 TcV. In particular, we predict the total pp cross section at the LHC to be 

116±4mb. 

Our results are summarized in Sect. 5, where a brief discussion is also presented. 
2. Regge analysis of tttt scattering (s 1 / 2 > 1.4 GeV) 



We normalize scattering amplitudes to 
4tt 2 

a AB = -TYjjj 2 T\ 1uiF a+b^a+b(s,0); X(a,b,c) = a 2 + b 2 + c? - 2ab - 2ac - 2bc. 

a ab is the total A + B cross section; for NN (pp, pp) and ttN scattering, we understand that the cross 
sections are spin-averaged. According to Regge theory, the imaginary part of a scattering amplitude with 
fixed isospin in the t channel, lmF A I ^ B ^ A+B (s,t), factorizes 3 as a product: for each Regge pole, R, we can 
write 

Imif;^ A+s (M) f { A R) (t)f { B R) (t){s/s) an{t) - (la) 

t fixed 

Here s is a constant, usually taken to be 1 GeV 2 ; we will do so here. A similar formula holds for the real 
parts: 

Reif;UA +S (M) Re£(R) f A R) (t)f B R) (t)(s/sr»W. (lb) 

t fixed 

£(i?), with Im£(i?) = 1, is known as the signature factor; for the Pomeron (P), P' and rho Regge poles one 
has 

M , „ 1 + cos7tq!b . „ _ „ 1 — cos7ra a , , 

Re£(i? = : -, R = P, P'; Re^(p) = — : p -. lc) 

sm 7ra/f sin na p 

The residue functions (t) depend on the quantum numbers of the Regge pole exchanged, on the particles 
that couple to it and, if we had external currents, also on their virtuality; but the power (s/s) aR ^ is universal 
and depends only on the Regge pole exchanged in channel t. The exponent ctn(t) is the Regge trajectory 
associated to the quantum numbers in channel t. For the Pomeron, which is rather flat, we will take it linear; 
for the rho, a more precise quadratic formula may be used. We thus write, for small t, 

a P (i) ( ~ o ct P (0) + a' P t, a p {t) ~ a p (0) + a' p t + \ol' p t 2 . (2) 
For the p and Pomeron pole, fits to high energy 7riV and NN processes give 

..' _ n nnri„u- 2 _ no r<„\r-i. 

(3) 



ap(0) =0.52 ±0.02, a' p = 0.90 GeV , a" = -0.3 GeV" 
a P (0)=l, a'p = 0.2 ± 0.1 GeV -2 . 



The Regge parameters taken here are based those in the global fit of Rarita et al.,1 13 ! which are still the 
best available as there are few modern data for the slopes in the relevant energy range. There are a few 
differences, however. For a p (0), we take the value 0.52 ± 0.02, instead of 0.58. This is more consistent with 
determinations based on deep inelastic scattering (see e.g. the paper of Adel et al., in ref. 12). Moreover, 
for Q.p(i) we use a quadratic formula that agrees with the average slope of ref. 13 for small, negative t, and 
which fulfills the condition a p (M 2 ) = 1. Finally, for a' P , Rarita et al. give 0.11, Froggatt and Petersen^ 0.3 
and the shrinking of the diffraction peak at the Tevatron suggests 0.26. Our choice here encompasses these 
three values. These are minor improvements as, in fact, for our fits in the present paper we only need the 
values of the q;_r(0); the slopes only intervene in sum rules. 

Let us now turn to the functions fi (t) . With respect to them we have two quite separate questions. 
First of all, we have the question of their normalization, that is to say, the values fi(0). These can be obtained 



3 In potential theory factorization can be proved rigorously; in relativistic theory, it follows from extended unitarity 
or, in QCD, from the DGLAP formalism. 112 ' 
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with little ambiguity and small errors by fitting experimental NN, irN and tttt total cross section data; we 
will do precisely that below. A different matter is the dependence of the fi(t) on t, i.e., the ratios fi(t)/fi(0), 
which is important in particular for Roy equations or sum rules like the ones at the end of the present 
Section. These are obtained from fits to the slopes of NN, ttN differential cross sections. Unfortunately, 
these fits are not unique, because both the background and the functional forms assumed for the fi(t) have a 
nonnegligible influence on the results, and because for the differential cross sections also the real part of the 
scattering amplitudes intervene. Moreover, the parameters of these fits were obtained before QCD emerged 
as the theory of strong interactions; these fits were extended to large values of t where, as we now know, 
Regge theory must fail and one has instead the Brodsky-Farrar behaviour. I 14 l They are thus forced fits. 

The situation, however, is not hopeless; the difference between the numerical results of various fits 
is small, for small values of \t\. For example, the numerical difference for the ratios fp(t)/ fp(0) between 
refs. 9, 16 is below the 10% level for It] 1 ' 2 < 0.4 GeV, which covers the values of t in which we are interested 
here. In the present paper we have chosen the t dependence of ref. 13, which was obtained in a detailed fit 
to many data. 

Before writing explicit formulas for the various processes (NN, irN, tttt) we have to decide in which 
variable we assume Regge behaviour to hold, which is important for us since we are going down to rather low 
energies. In (1) we have taken the cm. energy squared, s = (pi + P2) 2 , with pi the momenta of the incoming 
particles. Other possibilities are the s-u crossing symmetric variable, v — 2pi ■ p 2 , and E^ in , so we could 
assume behaviours like v ap or E 2 " n p instead of s ap , etc. We have, in our fits, tried all three possibilities; the 
fits using s, as in (1), all have substantially better x 2 /d.o.f. than those using v = 2pi or E^ in . Therefore, 
we stick to Regge behaviour in the variable s, as in (1). 

Regge formulas for tttt, ttN and NN scattering. We start with nn scattering. For exchange of isospin I t = 
in the t channel, containing the Pomeron and P' pole (the second associated with the /2(1270) resonance) 
we have 

Imi^=°)(M) s5L P(s,t) + P'(s,t), 

t fixed 

P(M)=frap(t)l±^^/r W , 

2 (4a) 

ap'(0)[l + ap'(0)J 
6 = (2.4±0.2)GcV~ 2 . 

Here (3 P = [fP] 2 , (3 P , = [fi^] 2 . 

The expression (4a) is like its counterpart in ref. 13, except for the P' pole parameters. In fact, 
the subleading contribution of the P' pole, that is necessary at the lowest energy range, is added somewhat 
empirically; its parameters are not well known, and we simply assume the corresponding trajectory to be 
degenerate with the one of the rho, as is suggested by a number of theoretical developments (in particular the 
QCD theory of Regge trajectories' 12 !), and as is done in ref. 6: ap'(i) = a p (t). In ref. 13, a larger value (0.7 
instead of 0.52) was given for the intercept of the P' pole and a smaller number was taken for its residue, but 
more modern determinations, as well as our fits, substantiate our choice; see Sect. 4, where we will present 
a global fit to data leaving, in particular, op'(0) as a free parameter. The result for it, ap>(0) = 0.55±0.03, 
is in perfect agreement with other modern determinations and with the degeneracy assumption. 

It should perhaps also be remarked that Eq. (4a), in what respects the Pomeron, is of limited validity 
(up to 10 - 15 GeV) since, at higher energies, total cross sections are known to rise. A modification of P(s, t) 
in (4a) that will make the parametrization valid up to multi-TeV energies will be given in Sect. 4. 

For I t = 1, we also take the parametrization of ref. 13. We write 

TmFg* =1 \a,t) t); 



t fixed 



1 1 + M*) 



p(s,t)=(3 p (1.5 + l)c w -1.5 " ' (s/s 

1 + a>p(Q) 



(4b) 



b is as before and (3 p = [/i- p '] 2 . The universal value of the slope of the diffractive factor, e bt , for all three 
trajectories rho, P and P', is what was found in ref. 13 from fit to actual NN and ttN data; it can nowadays 
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FIGURE 1. The average cross sec- 
tion |[2cr lr o, r+ + tr x o x o], which is 
pure It = 0. Continuous lines, for 
s 1/2 > 1.4 GeV: Regge formula. The 
lines cover the errors in the values of 
the Regge residues. Continuous lines, 
up to s 1//2 = 1.4 GeV: experimental 
cross section (from the fits in ref. 11; 
actually, with a slightly improved D2 
wave). The dotted and dashed lines 
are representative of the experimental 
errors in the cross section. 



be understood physically as a consequence of the universality of the Regge mechanism in QCD. We note 
that Froggatt and Petersen, I 6 ' who fit 7r + 7r~ data, find a value for b similar to ours for the Pomeron, but 
somewhat different ones for rho and P' . This last fact is not very meaningful as, in the fits to 7r + 7r~, the 
p, P' Regge poles are subleading and easily hidden by the Pomeron. We also remark that, in ref. 11, we 
had added a small background to lm.F^~ ^ to join smoothly the asymptotic formulas to the experimental 
cross section at s 1 / 2 <~ 1.4 GeV. With the value of the parameter (3 p found in the present Section, such 
background is unnecessary. 

For nir scattering we have to add an amplitude for exchange of isospin 2, corresponding to double 
rho exchange, which we do by writing 

\mF^= 2 \s,t) * R 2 {s,t) = p 2 e bt {s/s) a ^ +a ^- 1 . (5) 



We will discuss this quantity R2(s,t) later on; in particular, we will determine the quantity #2, which is 
small. We will start by putting f3 2 — and correct for this afterwards. 

The important parameters are (3p, (3p>, /3 p . We can obtain them fitting NN (pp and pp) and irN 
cross sections (including the forward differential cross section for the charge exchange reaction ir~p 
from 7T7T cross sections or from a global fit to the two sets. We write 



7r°n) 



da(ir p 



'pp 



TT°n) 



'pp 



Here f N/% 



with 



JN 



At 



s largo A 1 / 2 (s, nip, m p ) 
4tt 2 

' s largo A 1 / 2 (s, m 2 , m 2 ) 
2 1 — cos na 

JN/tt 



If 2 

2 JN/tt 



P(s,0) + P'(s,0) 
1 



V6 

^3 



4=0 s large 

( p ) / f ( p ) anc j we have defined 



sin 2 -kol l 



A(s,m 2 ,m 2 ) 



P(s,0) + P'(s,0) T^( s .°) 



l/5M)r 



(6a) 



p(s,t)=p p N ^ [(1.5 + l)e w -1.5 



1 + a p (t) 
l + a p (0) 



(s/S) 



<x P (t) 



o(Ntc) _ [ ,(P) Ap) 1 f (p) AP)] a 



(6b) 



(6c) 



In Eq. (6a) we have put the same values of /jv/w for Pomeron and P'. In Sect. 4 we will discuss fits allowing 
for different fffi n , /^J; their central values will be somewhat different, but the improvement in the x 2 /d.o.f. 
obtained by so doing is not significative. 

Fits. We will not fit data for scattering off neutrons which would not improve the precision while, because 
the neutrons are necessarily bound, they could distort the fits. We will also not include the difference of 
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cross sections a Pp — a pp in the fits, as this would involve the contribution of at least three Regge poles (ui, 
(j) and 7r) which do not contribute to irir. One could include the reaction pp — > fin, which only involves 
exchange of the rho, but the data for it are few and with (comparatively) large errors, so it would add 
little to the analysis. For the charge-exchange reaction, n~p — > 7r°n, only data in the forward direction are 
included. This reaction is interesting in that, although it has much larger errors than the others, it receives 
contribution from the real part of the corresponding Regge pole, so it represents a completely independent 
test of the Regge formulas. 

Before going on to the actual fits, a few words have to be said on the energy regions in which one 
may expect Regge behaviour (and, in particular, factorization) to hold. Generally speaking, we expect this 
to occur when one is past the region of elastic resonances and one also has E 2 in ^> A 2 (A ~ 0.4 GeV is the 
QCD parameter), which means for -Ekin ~ 1 GeV; but the precise details vary for different reactions. Thus, 
for pp, pp scattering there are no resonances and hence Regge behaviour is expected to occur precociously: 
here we will actually fit from Skin = 0.98 GeV. 

For tttt scattering it is difficult to tell when exactly one may use Regge formulas since data, par- 
ticularly for 7r~7r~, arc not very good. For the cross section cj( It=0 '> = ^[2a v o„+ + a„o„o], Eq. (4) provides 
a good representation for energies as low as Ekin = 1 GeV, as shown in Fig. 1; but, when resonances are 
more important, Regge behaviour is a good approximation only at slightly higher energies. Another matter 
is that, at low energies (s 1 / 2 ~ 1.5 GeV) the tttt data are of poor quality. Because of this, we will consider 
two extreme possibilities for actual fits. The first, that we will call no-cut, consists in including all tttt data 
for E kin > 1.1 GeV (s 1 / 2 > 1.38 GeV). The second possibility, that we call cut, consists in cutting out all 
data for energies below s 1 / 2 = 2 GeV. The difference in results between the two fits will be an indication of 
the systematic errors in our calculation. 
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FIGURE 2. The total cross sections <J n ± p , \(p PP + a vv) an< i ^( a K+ P + °K-p)- Black dots, trian- 
gles and squares: experimental points. Continuous lines: Regge formulas, with parameters as in our 
best fit. For \(<Jp p + cr PP ) and \{u K + p + cr K - p ), the three lines cover the errors in the values of the 
Regge residues. For ttN the theoretical error is of the order of that for ^(<Jp P + a pp ). Note that the 
thick line in the low energy experimental cross sections for nN is merely due to the accumulation of 
closely spaced data. 
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For ttN the formulas (6) fit well data down to ~ 1.3 GeV, but, for the sum, cr„.+ p + cr„- p , one 
can go to -Ekin ~ 1 GeV. For the difference, cr n + p — cr n - p and for the charge-exchange reaction tt~ p — > ir°n, 
resonances somewhat spoil local agreement, but Eq. (6) provides a good average representation even down 
to 1 GeV, as has been known for a long time (see, e.g., ref. 15) and as can be seen in the lower energy region 
in our fit to ir + p data in Fig. 2. We will here start from £^in = 1-08 GeV. 

Another question is how high one goes in energy. In the present Section we fit experimental data 
for cm. kinetic energies Skin ^ 16.5 GeV: this is what is required for applications to tttt Roy equations, 
dispersion relations and sum rules, since here the importance of the very high energy region is negligible. 
Nevertheless, and as stated before, parametrizations and fits valid up to multi TeV energies will be given in 
Sect. 4. 

The data on ir~p — > ir°n are from the compilation in rcf. 15. For NN and ttN we will take the 
data from the COMPAS Group compilations, as given in the Particle Data Tables. I 16 ' For those data where 
systematic errors are not given, we have included a common systematic error of 0.5% for pp, 1% for pp and 
1.5% for irp, which are like the standard systematic errors in other data. Another possibility is to take a 
common systematic error of 1.5% for all data: the difference of the results with the two will indicate the 
systematic errors of our fit. We have also made a selection of NN, irN data, as follows. We take only data 
at energies at which there are results for both pp and pp or ir + p and tt~p; and, when there are, at a given 
energy, data from various experiments, we have taken only the most recent. This is designed to thin out 
the data to a number comparable in order of magnitude to that of tttt, so that nn data have a nonnegligiblc 
weight in the joint fits. For tttt scattering we have taken the errors as given by the various experimental 
groups except for those of Abramowicz ct al.,^ who only give statistical errors, much smaller than those of 
the other groups, and for which we have added a common systematic error of 1.5 mb to all points; even with 
this, the error, though comparable, is smaller than what other groups find. 

We could fit separately the NN, ttN data and the tttt data of ref. 7, or make a global fit. The results 
of these fits, in which we have put 02 = 0, and fixed a p (0) = 0.52, are given in Table I, where the errors 
correspond to one standard deviation. The best values are average values, with errors enlarged to overlap 
other results. A graphical representation of this best fit may be seen, compared with experimental NN, ttN 
cross sections in Fig. 2, and for tttt data, in Fig. 3. We note that, in Fig. 3 for tttt, we have used the values 
of p and 02 from Eq. (11) below. 





NN, ttN [enlarged error (a) ] 


Only tttt [cut (b) (c) ] 


NN, ttN, tttt [cut (b) ] 


Best values 




1.405 ± 0.005 [1.411 ± 0.004] 




1.404 ± 0.005 [1.407 ± 0.005] 


1.406 ± 0.007 




0.366 ± 0.009 [0.367 ± 0.010] 




0.366 ± 0.010 [0.367 ± 0.009] 


0.366 ±0.010 


Pp 




1.30 ±0.13 [0.59 ±0.27] 


1.32 ±0.13 [0.59 ±0.25] 


1.0±0.3 W 


Pp 


2.55 ±0.01 [2.53 ±0.01] 


2.50 ± 0.08 [2.55, fix] 


2.56 ±0.01 [2.56 ±0.01] 


2.56 ±0.03 


f3 P , 


1.04 ±0.02 [1.09 ±0.02] 


1.46 ±0.17 [1.04, fix] 


1.04 ±0.02 [1.04 ±0.02] 


1.05 ±0.05 


x 2 


303 r 252 "I 


109 r 45 1 


415 r 348 ] 




d.o.f. 


229 - 4 L229 — 4 J 


58-3 L 39 — 1 J 


288 - 5 1.268 -5J 



< - a ' We here endow all ttN numbers with a minimum systematic error of 1.5%. ' By "cut" we mean that tttt data 
for s 1/2 < 2 GeV are removed from the fit. We here fix (3 P , P i as given by NN, ttN to avoid spureous minima. 
The error in this quantity will be improved using crossing sum rules; see Eq. (11) below. 

Table I 
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A few features of our results worth noting are the following. Firstly, the equality of f N j v and (3p, 
f3pi , for fits with and without tttt data is a very satisfactory test of factorization. Another interesting point 
is the stability and accuracy of the parameters /jv/ttj P p N ^\ Sip- The parameter [3p> is less well determined, 
and (3 p is not fixed with precision by fits to data alone; we will improve its accuracy in a moment using 
sum rules. Secondly, the matching between the low energy (s 1 / 2 < 1.42 GeV) results for cross sections from 
phase shift analyses, and the high energy (s 1 / 2 > 1.42 GeV) Regge representations is excellent for 7r°7r~, 
tt~tt~ and cr( /t= °). It is less good for tt + tt~ , where matching occurs only at the 1.5 a level, no doubt due to 
the coinciding tails of the /2(1270) and /o(1370) resonances. And, thirdly, the fact that, for NN and ttN 
the x 2 /d.o.f. is somewhat larger than unity is due to the following effects. First, we use only two poles for 
vacuum exchange, and one for charge exchange: we are thus missing the contributions of other poles, likely 
small, but not negligible at the lower energy range. Secondly, at the very low energy range, the experimental 
cross sections oscillate a little around the Regge formulas. Finally, and probably the most important effect, 
we have that, to cover well the upper part of the energy range, we need more sophisticated expressions: see 
Sect. 4. 

Besides this, we have a few technical points to make in connection with the fits including tttt data. As 
is clear from Fig. 3, the low energy (s 1 / 2 < 2.5 GeV) results for tt~tt~ cross sections of various experiments 
are quite incompatible with one another, which is the reason for the large x 2 /d.o.f. in no-cut fits. There is 
certainly a bias in the experimental tt~tt~ cross sections of Biswas et al., and Robertson, Walker and Davis, ^ 
in the lower energy range. This is probably due to incorrect treatment of final state interactions, that, at 
these lower energies, are influenced by the A33 and other resonances. At higher energies the influence of 
this resonance seems to become negligible as, indeed, the n~ir~ cross sections found by Robertson, Walker 
and Davis overlap those of Abramowicz et al.' 7 ' and both tend to the tt + tt~ one, as Regge theory and the 
Pomeranchuk theorem imply. We consider that this problem is solved by considering our two types of fits, 
cut or no-cut, for tttt scattering. 

We next discuss the isospin 2 exchange piece, i?2(s,f). We have three methods to get the quantity 
/?2- First, we fix the values of (3p and f}p> to their best values, as given in Table I, and fit the tttt data using 
Eqs. (4), (5). Note that one cannot leave the parameters (3p, [3p> free in these fits because one would get 
spureous minima, since the data are not precise enough. We find (3 P = 1.07 and a very small /3 2 ~ — 2 x 10~ 8 . 
Alternatively, we could obtain fa by fitting cr n o v o — cr 7r o 7r + at s 1 / 2 = 1.42 GeV, as was done in ref. 11. This 
gives 02 — 0.55 ± 0.2. Finally, we can use the first crossing sum rule in the Appendix to ref. 11 (identical to 
(B.7) in ACGL), which would give a [3% compatible with zero. We take as a compromise the number 

/3 2 = 0.2±0.2. (7) 

However, we should note that the t dependence of i?2(s,i) is little more than guesswork. 

Sum rules. We now say a few words on the sum rules discussed in ref. 11. Because these sum rules were 
verified with Regge expressions slightly different from what we have now found, one may wonder what 
happens to them. Since the formulas in (4), (5), with parameters as in Table I, agree with those of ref. 11 
within ^ 2 a, and the decrease of (3p is (partially) compensated by the increase in [ip> , it can be expected 
that the various sum rules would still be satisfied within errors, as indeed it happens. Our numbers here 
leave the agreement of the Olsson sum rule and the value of the P wave scattering length and effective range 
still within 1 a. We have already discussed the first crossing sum rule in the Appendix in ref. 11 in connection 
with /3 2 , so we turn to the second crossing sum rule. It reads, 



r ds ^' J - 1) (M^)-M(H M)= ,« d /^[ g -2rJ Im 

Jiml s J 4ml S 2 (s-Amiy 



(8) 



The interest of this sum rule lies in that its high energy (s 1 / 2 > 1.42 GeV) is dominated by p(s, t), while the 
low energy piece (s 1 / 2 < 1.42 GeV) is such that the contributions of the S waves cancel, so it is dominated by 
the P wave, which is very well known. Thus, it provides an independent, reliable way of fixing the parameter 
[3 p . We find (8) satisfied provided one has 

P P = 0.82 ±0.12. (9) 
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FIGURE 3. Total cross sections o n a n -, o^-^- and a^+ n -. Black dots, triangles and squares: ex- 
perimental points from ref. 7. The stars at 1.38 and 1.42 GeV (PY) are from the phase shift analysis 
of experimental data given in ref. 11, slightly improved for the D2 wave. Continuous lines, from 1.42 
GeV (PY): Regge formula, with parameters as in our best fit (the three lines per fit cover the error 
in the theoretical values of the Regge residues). Dashed lines, above 2 GeV: the cross sections follow- 
ing from ACGL;' 1 ' the grey band covers their error band. Below 2 GeV, the dotted line corresponds 
to the n + n~ cross section from the Cern-Munich analysis; cf. Fig. 7 in the paper of Hyams et al.' 8 ' 
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[Our best values] 


[PY] 


[ACGL] . 


Pp 


0.94 ±0.14 


0.84 ±0.10 


1.48 ±0.25 


Pp 


2.56 ±0.03 


3.0 ±0.3 


1.0 ±0.6 


Pp> 


1.05 ±0.05 


0.72 ±0.07 


2.22 ±0.38 


P2 


0.2 ±0.2 


0.55 ±0.20 






Since this is compatible with the independent determinations in Table I, we may include fulfillment of (8) 
in the fits. If we do so for the fit with cut tttt data, we get the value 

P p = 0.78 ±0.11. (10a) 

If we include (8) in the fit with all 7T7t data (no-cut) we find instead 

P P = 1.07 ±0.09. (10b) 

Combining (10a, b) we can then take 

P P = 0.94 ± 0.10 (Stat.) ± 0.10 (Syst.). (10c) 

-Best values. We can now present our best values, and compare them with the values given in ref. 11 (PY), 
obtained basically from those by Rarita et al.,' 13 ' or those of rcfs. 1, 5 (ACGL): 



(11a) 



Besides these, we have also 

f N/ „ = 1.408 ±0.005, P p Nn) = 0.366 ± 0.010. (lib) 

Our present results are compatible with those in refs. 6, 11, 13. We note, however, that our fits 
include much more information on the total cross sections than those in refs. 6, 13. The first only includes 
7r + 7r~ data while the more complete fit of Rarita et al.I 13 ' includes 24 total cross section data for NN (we 
have 34) and 28 for ttN (we have 87); the energy range we cover is also wider, by a factor 6 in the variable s. 
We also have 58 tttt data points (none in ref. 13). Of course, the situation is different for the t dependence 
of the residue functions /,(t) for which the fit of Rarita et alJ 13 ' cannot be really improved. 

The results in (11) may be compared with some theoretical models. The value f N /^ ~ 1.4 is 
similar to what one gets in the naive quark model' 17 ! with additive quark-quark cross sections, that gives 
/jv/tt = 3/2. [It is, however, not clear why the naive quark model works, as its mechanism is very different 
from the orthodox QCD one]. Likewise, the value of j3 p = 0.94 ± 0.14 is similar to what one has in the 
Veneziano model! 18 ! (^ ~ 0.95), but the relation (3 pi = that this model gives is not well satisfied. P p 
also agrees with the rho dominance model, in which one couples the rho universally to pions and nucleons 
according to 



with t = (?/2, (t the Pauli matrices, that gives P p = ^/| /jv/tt Pp N ^ — 0.84. 
3. ttK scattering 

The analysis of ttK scattering follows similar lines. For exchange of isospin zero we have 

Im i#T°> (a, t) s5L f Kh [P(s, t) + rP'(s, t)] ; 

t fixed 

/K/.=/r(o)//i p) (o). 



(12a) 
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P, P' are as above, and r is related to the branching ratio for the KK decay of the resonance 4 /2(1270), 
r <~ BR = (4.6 ± 0.5) x 10~ 2 . For isospin 1 exchange, 

t fixed (12b) 

p(s, t) is as before. To find the desired representations for the nK amplitude we have to determine the ratios 
Jk/tt, 9k/tt- For the first, this is done taking the /jv/tt from NN, nN scattering, as in the previous sections, 
and with the help of the even combination of cross sections for KN scattering: 

4tt 2 r i 

*k+ p + *k- Ps £ ge x y2 {s ^ ml) fN/J Kh [P(s, 0) + rP'(s, 0)j . (13) 

For gK/TT, unfortunately, we cannot use the charge exchange reaction K~p — ► K°n because there are two 
trajectories of comparable importance, p and that corresponding to a2(1320) exchange, that contribute; for 
a discussion, cf. for instance the text of Barger and Cline, ref. 12. The difference of cross sections K + p and 
K~p also contains extra contributions (u>, <f), . . .). 

For the KN cross sections we will take data in the region E kln > 1 GeV, and go up to E^ in = 10 GeV. 
At higher energies the logarithmic increase of the total cross section for K + p scattering is noticeable, and we 
would need more complicated Regge formulas (that we will give in Sect. 4); while, as occurs for the 7T7t case, 
the importance of the very high energy region is negligible in most applications to nK scattering. For nK 
scattering we thus expect the ensuing Regge expressions to be accurately valid for a corresponding energy 
range, say, for 1.7 GeV < s 1/2 < 11 GeV. 

The K ± p data we take also from the COMPAS Group compilations; see the Particle Data Tables. I 16 ' 
For those data where systematic errors are not given, we have included a common systematic error of 0.3 mb, 
as we did for the nN case. We take only data at energies at which there are results for both K + p and K~p. 
In the fits we use the very precise values of the parameters /at/^, fip obtained before, and we set r = 0, since 
it is very small and not very well known. In fact, in Sect. 4 we will make fits leaving r free; its value will 
turn out to satisfy \r\ < 0.1. We find 

f K/7r = 0.67 ± 0.01 [from K+p + K~p; X 2 /d.o.f . = 50/(43 - 1)] 

9K/TT =1.1 ±0.1. 

The results for (a K + p + a K - p )/2 are shown in Fig. 2. The value of g K / v is taken from the classical analysis of 
ref. 19, that takes into account the a2(1320) exchange. The value of /k/w is within 20% of its SU(3) value, 



A global fit valid up to multi TeV energies 



A simple parametrization of scattering amplitudes which fits data at energies s 1 / 2 > 12 GeV (with a 
X 2 /d.o.f. = 1.2 to 1.8, depending on the process) may be found in in refs. 20, 21. Here the Pomeron is 
allowed an intercept larger than unity, ap(0) <~ 1.095, and the intercept of the P' is given as ap'(0) = 0.66. 
This parametrization, that we will call "power Pomeron" parametrization, is purely phenomenological, as 
explicitly mentioned in refs. 20, 21. Only data with energy larger than ~ 10 GeV are used in the fits 
which, if extended to energies below 5 GeV, miss widely the data. These parametrizations also must fail 
at very large energies since they are incompatible with unitarity in that they violate the Froissart bound. 

4 Since the P' pole couples so weakly to kaons, one may wonder on the importance of other Regge poles for the 
subleading contribution to kaon scattering. For KK scattering, the Regge pole associated with the /2(1525) 
resonance gives a substantial contribution; but, for KN or nK scattering, this trajectory contributes very little 
since it is almost uncoupled to pions and nucleons and its intercept is small, a/ 2 (i B2 B) ^ —0.3. For KN and nK, 
the amplitude for exchange of zero isospin is almost pure Pomeron. 
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As a matter of fact, in rcf. 22 it is remarked the inadequacy of such parametrization, and a parametriza- 
tion verifying the Froissart bound (i.e., with a term in (Const.) x log 2 s/so + Const.) is substituted in 
place of the "power Pomeron." This improves substantially the x 2 /d.o.f. of the fit, and gives an intercept 
ctP'(O) = 0.54 ± 0.02, perfectly compatible with our choice 0.52 ± 0.02. The corresponding parametrization 
holds down to s 1 / 2 = 5 GcV. 

It is possible to write a parametrization, similar to that of ref. 22, obtained by a modification of 
the Pomeron in Eq. (4a), that fits data for kinetic energies from 1 GeV to the multi-TeV region and which, 
moreover, is compatible with unitarity by adding a slightly more complicated logarithmic term. We do this 
as follows: we note that it is possible to improve the Froissart bound to a bound of the form! 23 ! 

cr tot < a log 2 5 , (15) 

si log ' s/s 2 

which is maximal in the sense that one cannot increase the power of the log in the denominator to more than 
7 / 2 . For the bound for tttt scattering, one can evaluate the constants a, si, s 2 in terms of the pion mass and 
low energy parameters for the D wave, with a = 7r/4m 2 ~ 15 mb 2 , s\ = m 2 if we assume the cross section 
to be mostly inelastic. What this suggests is that we add a term like (15) to the Pomeron given in (4a), but 
leaving a, s\, s 2 as free parameters. Thus we replace, 

P(s,t) =(3 P a P (t) 1 + a 2 p{t) e bt ( S /§r^ P F (s,t), 

P F { S ,t)=k P + A\oi * )a P (t) 1 + ap{t K bt (s/sr^. 

{ si log 1 s/s 2 J 1 

This replacement should also be made in Eqs. (6), (12) and (13). The logarithmic term has an appealing 
physical interpretation as the contribution of the Regge cuts which, as Mandelstam showed long ago,' 24 ' 
should accompany the J'omeron. The parameter (3p that we used before is to be viewed as an effective 
parameter, the sum of f3p and the average value, for low energy (s 1 / 2 < 15 GeV), of the logarithmic piece 
in (16). 

With (16) we fit data for ir^p, K + p + K~~p, tttt and pp + pp cross sections 5 up to the highest energies 
attained experimentally, 30 TeV in cosmic ray experiments. I 25 ' 

Because we have so many experimental data, covering such a wide energy range, we may fit all 
hadronic data (i.e., including NN, ttN, KN and tttt data) leaving all parameters free; in particular, this will 

test the quality of the assumption of degenerate rho and f 2 trajectories, the equality of fff n , /^J, and the 
smallness of the parameter r in Eq. (13). We only fix a p (0) to the number given by deep inelastic scattering, 
0.52 ± 0.02, include the sum rule (8), and find 

1.350 ±0.008; /<£> = 1.67 ± 0.07; f K/ „ = 0.74 ± 0.01; a P ,(0) = 0.61 ± 0.04, 
P P =2.33 ±0.09; (3 P , = 1.05 ± 0.10; (3 { p Nw) = 0.385 ± 0.009; P = 0.94 ± 0.14; r = 0.11 ±0.08, 
A =0.022 ±0.002; si = (1.2 ± 0.7) x 10~ 4 GeV 2 ; s 2 = (0.33±g;|) x 10~ 7 GeV 2 , 
X 2 /(d.o.f.) = 339/(358 - 12) ~ 0.98. 

(17) 

The value of f3 p given here is that found before, Eq. (10c); since there are no hit data at very high 
energy, the value of this quantity essentially decouples from the very high energy analysis. 6 

What is interesting of (17) is that the value of a P '(0) is compatible with what one finds from 
degeneracy, a P '(0) = a p (0) = 0.52±0.02, and that f N L and f N ^ are not far from each other, as required by 



Above 30 GeV we approximate op p — a pp — (66.7mb)(s/s) -0 ' 55 , where this difference comes from the phenomeno- 
logical fit of ref. 16, since we do not have data at coinciding energies. For tttt only data above 2 GeV are included 
in these fits. 

If we had fitted also f3 p , including the sum rule (8), its value would depend on whether we had included all nn 
data above 1.4 GeV (in which case we would have got 1.05 ± 0.009) or only data for s 1 ^ 2 > 2 GeV, which gives 
0.80 ± 0.11: essentially the same numbers as in the fits in Sect. 2, Eq. (10a,b). 
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FIGURE 4. The total cross sections a 7T ± p , \{cr K + v + a K - p ) and \{<Jp p + u pp ) up to 30 — 60 GeV 
(upper graph) and \{crpp + cr pp ) up to 30 TeV (lower graph). Black dots, triangles and squares: 
experimental points. For energies above 30 GeV, we have given the experimental values of \{crp P + 
a pp ) as if they equaled a Pp or a vp . Continuous lines: Regge formulas, with parameters as in (18). 
In the lower figure we have given the error bands for \(p pv + a pp ) that follow from (18). 



(strong) factorization. In fact, this had alredy been noticed in rcf. 22: in a fit with a formula compatible with 

theory (the Froissart bound), the resuts respect other theoretical constraints automatically. The problem 

with the fit in (17) is that there is, unfortunately, a very strong correlation among 0p, 0p>, ap<(0), si and 

S2 and, if we leave all of them free as we did in getting (17), there exist a large number of equally significant 

minima: the parameters are not well determined. In fact, si, s 2 , 0p< and ap'(O) can one mock the effects 

of each other. In particular, a set of fits with quality essentially unchanged may be obtained by varying 

(p) (p r ) 

simultaneously Si and s 2 . In view of this, we require = f N /^ and , to fix the parameters, choose 

Si = 0.01 GeV 2 and repeat the fit with all other parameters free. We find what we consider our best result: 
f N/n =1.380 ±0.004 f Kh = 0.717 ±0.005; a P -(0) = 0.55 ± 0.03; r = ± 0.013, 
P =2.31 ±0.05; P , = 1.39 ± 0.14; p N ^ = 0.377 ± 0.009; p = 0.94 ± 0.14, 

(18) 

A =0.033 ±0.001; si = 0.001 GeV [fix.]; s 2 = 0.13 ± 0.05 GeV ; 
X 2 /(d.o.f.) = 372/(358 - 10) ~ 1.066. 

We note that, although the x 2 /d.o.f. is slightly worse than that in (17), we consider the fit in (18) to be 
more satisfactory physically. The values of the parameters s\, s 2 in (17) were too small for comfort, and one 
should not force too good a fit at the expense of physical considerations (like factorization or degeneracy), 
particularly since we are fitting with formulas that, at the lowest energies, should be corrected by including 
other Regge poles (or cuts). Eq. (18) has the nice properties that degeneracy is verified, up to errors, that 
Jk/tt agrees better with its SU(3) value, and the value of 0p* agrees, also within errors, with the prediction 
of the Veneziano model, 0p> = \0 P — 1.4 ± 0.3. 

At the lower energies (below 15 GeV) (16) plus (17) or (18) overlap with the previous fits, using 
(4a) for the Pomeron and P', for vacuum exchange. In fact, for Kp or nN, the corresponding curves could 
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not be distinguished from those obtained using (4a) in Fig. 2; see Fig. 4. For pp + pp, the result of the fits 
with the two types of formulas, (4a) and (18) are depicted in Fig. 4, where the error bars corresponding to 
(18) are also shown. 

5. Summary, and a short discussion 

The Regge parameters that ACGLW and, following them, the authors in refs. 2, 3, 4, 17 assume not only 
are unorthodox but, as we have shown, incompatible with experiment. As our Fig. 3 clearly demonstrates, 
the claimed large errors in ACGL are not large enough to cover the experimental data. 

ACGL get these quaint Regge parameters by considering sum rules like (8) that link the Regge 
contributions, which they assume to hold only for s 1 / 2 > 2GeV, with the corresponding low energy (s 1 / 2 < 
2 GeV) pieces. Unfortunately, the intermediate intermediate energy (1.4 GeV < s 1 ' 2 < 2 GeV) that ACGL, 
again here followed by the authors in refs. 2, 3, 4, take for the SO, P, DO and F phases come basically from the 
experimental analysis of the Cern-Munich group, whose tt + tt~ cross section is more and more incompatible, 
as s 1 / 2 nears 2 GeV -in fact, as soon as inelasticity becomes important- with the values found by all other 
experiments: ^ see our Fig. 3. [The interested reader may consult ref. 9 for the detailed discussion of this 
and other related issues]. It is thus not surprising that Pennington ^ and Ananthanarayan et al.,M who fix 
their Regge parameters by balancing them above 2 GeV with phase shifts below 2 GeV, get totally incorrect 
Regge amplitudes. And, given these facts, it also follows that the low energy results of references 2, 3, 4, 
which borrow their input at energies s 1 / 2 > 1.4 GeV from ACGL, should be taken with great caution. 

Unlike the results of phase shift analyses, the Regge formulas in Eqs. (4), (16) with the parameters 
as the "Best values" in (11) or (18), and which we summarize in Table II, give a consistent representation 
for the imaginary part of all the tttt scattering amplitudes, a representation which can be trusted, within the 
given errors, for s 1 ' 2 > 1.4 GeV, provided |i|V2 < 0.4 GeV. In fact, one has better than that: our Regge 
formulas give a good representation of those processes in pion-pion scattering where resonances are absent, 
or are not important, down to lower energies, just as it happens in NN or ttN scattering. This occurs, 
in particular, for ir°ir + and n~ir~ , for which the Regge formulas reproduce the experimental data down to 
s 1 / 2 ~ 1.1 GeV. However, by the very nature of things, we are likely to have uncertainties of the order of 
15% in the region 1.4 GeV < s 1 / 2 < 1.8 GeV when exchange of isospin 1 is important, because the Regge 
formula probably represents data only in the average there, as occurs for ttN scattering. Finally, and using 
Eqs. 6, 13 and the formulas in the last column in Table II, we can fit NN, irN and KN up to multi TeV 
energies, and predict tttt and ttK cross sections there. 

When performing calculations of tttt scattering in which the lower energy region is dominant (such 
as Roy equations, dispersion relations or sum rules) it is irrelevant, within our errors, which form one uses 
for the Pomeron, (4a), (17) or (18). The last has better overall fit, and (probably) a more realistic value for 
(3p>; although the first is to be preferred in that it is simpler and fits slightly better the low energy data. 
The safest procedure is to use both fits, and consider their difference as a measure of the influence of the 
parametrization in the results. We should, however, emphasize that the parameters in the fits are strongly 
correlated and, even when they are similar, one cannot mix parameters from the various columns in Tabic II; 
each fit stands on its own. 

One may also wonder what happens for values of the momentum transfer larger than It] 1 / 2 ~ 
0.4 GeV. On general grounds, one expects Regge theory to work when s ^> A 2 , s 3> \t\ and in fact, as 
already mentioned, Regge representations for NN or 7riV become unreliable at large \t\. For example, the 
parametrizations of Rarita et alJ 13 ' and ref. 6 for f^ p \t) differ completely from one another already at 
—t = 0.23 GeV 2 , where the first changes sign. There is unfortunately no sure way out of this problem 
(which is further discussed in the second paper in ref. 23), and one has to admit that, for s 1 / 2 > 1.4 GeV 
and values of the momentum transfer \t\ > 0.15 GeV 2 , there is no reliable information on the pion-pion 
scattering amplitude -which, in particular, is an unavoidable cause of uncertainty for Roy equation analyses 
that require information for values of \t\ as large as 0.5 GeV 2 . 
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1 GcV < E kin < 15 GeV 
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0.11 ±0.08 


± 0.013 
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0.52 ± 0.02 [fix] 


0.52 ± 0.02 [fix] 


0.52 ± 0.02 [fix] 
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0.52 ± 0.02 [fix] 


0.61 ±0.04 


0.55 ± 0.03 
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2.31 ±0.05 
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0.033 ± 0.001 
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(1.2 ±0.7) x 10~ 4 GeV 2 


= 0.01 GeV 2 


S2 




(0.33±°;®) x 10" 7 GeV 2 


0.13 ±0.05 GeV 2 


X 2 /d.o.f. 




0.98 


1.07 



Table II 
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